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The contributions to the electron-electron-photon vertex of two-loop electroweak corrections are 
calculated. The relative correction to the parity-violating asymmetry of Mpller scattering for the 
case of 11 GeV electron scattered off the electron at rest is found to be about —0.0034 and should 
be taken into account at futnre experiment MOLLER at JLab. 


I. INTRODUCTION 

Nowadays high energy physics faces difficulties. The energies reached at existing facilities are almost at the limit of 
their technical possibilities and reasonable cost of the projects. Besides whole set of experimental information is in a 
confident agreement with Standard Model (SM) predictions, there are large number of indications that at this scale 
there is a possibility for ” new physics” to manifest itself. This is mostly the content of scientific programs for existing 
and future accelerators. But most probable option is that within accessible energy range the traces of new physics 
can be discovered as a small deviations from SM predictions. Revelation of new physical phenomena is only possible 
by comparing of detailed experimental result with model predictions. The aim of the present paper is to continue the 
elaboration of precise description of one of the most prominent process — Mpller scattering [l|. 

This process has wide and active interest from both experimental and theoretical sides for several reasons. It has 
allowed the high-precision determination of the electron-beam polarization at SLG i, SLAG i a, JLab 0 and 
MIT-Bates 0 (and as a future prospects — the ILG [ 3 ] and CLIC 0 ). The polarized Mpller scattering can be an 
excellent tool in measuring parity-violating weak interaction asymmetries 0 . The first observation of Parity Violation 
(PV) in the Mpller scattering was made by the E-158 experiment at SLAG [lol - [T3| . which studied scattering of 45- to 
48-GeV polarized electrons on the unpolarized electrons of a hydrogen target. It results at = —t = 0.026 GeV^ for 
the observable parity-violating asymmetry Apy = (1.31 ± 0.14 (stat.) ± 0.10 (syst.)) x 10“^ [l^ which allowed one 
of the most important SM parameters — the sine of the Weinberg angle sin^vv — to be determined with the best 
accuracy at that moment. 

The MOLLER (Measurement Of a Lepton Lepton Electroweak Reaction) experiment planned at the Jefferson Lab 
aims to measure the parity-violating asymmetry in the scattering of 11 GeV longitudinally-polarized electrons from 
the atomic electrons in a liquid hydrogen target with a combined statistical and systematic uncertainty of 2 % [13- 
[I3- With such precision any inconsistency with the SM predictions will clearly signal the new physics. However, a 
comprehensive analysis of radiative corrections is needed before any conclusions can be made. Since MOLLER’s stated 
precision goal is significantly more ambitious than that of its predecessor E-158, theoretical input for this measurement 
must include not only a full treatment of one-loop (next-to-leading order, NLO) electroweak radiative corrections 
(EWC) but also two-loop corrections (next-to-next-leading order, NNLO). Although, two-loop corrections to the cross 
section may seem to be small, it is much harder to estimate their scale and behavior for such a complicated observable 
as the parity-violating asymmetry to be measured by the MOLLER experiment. 

This paper is the part of our attempt to perform this theoretical input. The significant efforts have been already 
dedicated to one-loop EWC. A short review of the references on that topic is given in , where we calculated 
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a full set of the one-loop EWC both numerically with no simplifications using computer algebra packages and by¬ 
hand in a compact form analytically free from nonphysical parameters. One way to find some indication of the size 
of higher-order (two-loop) contributions is to compare results that are expressed in terms of quantities related to 
different renormalization schemes. In we provided a tuned comparison between the result obtained with different 
renormalization conditions, first within one scheme then between two schemes. Our calculations in the on-shell and 
Constrained Differential Renormalization schemes show the difference of about 11%, which is comparable with the 
difference of 10% between MS [2l[ and the on-shell scheme [2^ . 

The two-loop EWC to the Born cross section MqMq) can be divided onto two classes: Q-part induced 
by quadratic one-loop amplitudes ~ ^ and T-part - the interference of Born and two-loop amplitudes 

^ 2Re(Alo-^^) (here index i in the amplitude Mi corresponds to the order of perturbation theory). The Im¬ 
part was calculated exactly in (using Feynman-t’Hooft gauge and the on-shell renormalization), where we show 
that the Q-part is much higher than the planned experimental uncertainty of MOLLER, i.e. the two-loop EWC is 
larger than was assumed in the past. The large size of the Q-part demands detailed and consistent treatment of T-part, 
but this formidable task will require several stages. Our first step was to calculate the gauge-invariant double boxes 
24|. Vext step was to calculate the two-loop gauge invariant set of boson self energies and vertices function diagrams 
25l| . In this paper we took into account two photon emission mechanism in soft photon approximation. It is important 
to calculate hard photons bremsstrahlung contribution in accordance with MOLLER detector parameters. This work 
was partially done for one-loop EWC in [2g. Then we consider in the EWC arising from the contribution of a wide 
class of the gauge-invariant Feynman amplitudes of the box type with one-loop insertions: fermion mass operators 
[or Fermion Self-Energies in Boxes], vertex functions [or Vertices in Boxes], and polarization of vacuum for bosons [or 
Boson Self-Energies in Boxes]. Also, in this paper one can find extended literature review and all necessary details 
and notations useful for understanding of present paper. So, finally, this paper logically related to manuscript [T7I| . 
where we do the next step — we calculate the insertions of two-loop vertices to vertices (VV), fermion self-energies to 
vertices (FSEV) and double vertices (DV). 

This theoretical program must be completed by testing our two-loop results with the results of computer algebra 
packages (like FeynArts, FormCalc etc). This testing was already done for one-loop calculations, to do this for two 
loops results is the important task for our group to be done in the next future. 

The paper is organized as follows. In Section |n] we consider the asymmetry in Born approximation and introduce 
the basic notations. In Section m we calculate two Feynman diagrams with extra W and Z boson subgraphs (VV). 
Section m is devoted to the diagrams with lepton mass operators insertions (FSEV). In Section |V] complex vertices 
are considered (DV). And in Section |VT] we give numerical estimation of total effect of these contributions. 


II. BASIC NOTATIONS 


We consider the process of electron-electron elastic scattering, i.e. Mpller process: 

e(pi, Ai) -b e(p 2 , A 2 ) e{p[,X[) + e(p 2 , A 2 ), 


( 1 ) 


where Ai ^2 (A^ 2 ) 9'''® f^e chiral states of initial (final) electrons and pi ^2 are 4-momenta of initial electrons and Pi 2 
are 4-momenta of final electrons. The first measurement of parity-violating (left-right) asymmetry 
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( 2 ) 


in Mpller scattering was made by E-158 experiment at SLAG [lOKl^- In lowest order of perturbation theory in frames 
of QED the matrix element squared which is summed over polarization states of electrons has the following form: 


= 8(47ra) 


2 s'* -b -b 


(3) 


We use the notation for the kinematic invariants neglecting of electron mass to: 

s = 2 pip 2 , ^ = —2pip'i, u = — 2 piP 2 , s + t + u = 0. 


(4) 


Thus here and further we neglect the terms of order /s) since in MOLLER experiment it is expected that beam 
energy is Abeam = H GeV, that is s = 2TOAbeam ~ 0.01124 GeV^. Within the Standard Model one has additional 
contribution in Born approximation with Z-boson exchange which gives rise to polarization asymmetry A^: 


aQ - ^ A ® 

4 * - r, 2 ^( 0 )-■ 
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Figure 1: Two-loop vertices to vertices (a, b), fermion self-energies to vertices (c,d), double vertices (e-h). Photons are denoted 
by wavy lines, massive bosons - by dashed lines, electrons and neutrinos - by solid lines. Notations on diagrams show the type 
of particle {Z, W, v) and 4-momenta of them. 


where c = cos0 is the cosine of scattering angle 9 = in the system of center-of-mass of electrons, mw is the 

FF-boson mass and a is the so-called ’’weak electron charge” 


a — 1 — 4s^. (6) 

Now let’s recall that Sw (cw) is the sine (cosine) of the Weinberg angle expressed in terms of the Z- and FF-boson 
masses according to the Standard Model rules: 


sw = \J^-c^, cw = 'mwlmz. (7) 

Thus, the factor a is just a ~ 0.109 and the asymmetry is therefore suppressed by both slrri^ and a. Even at Central 
Region (CR) of MOLLER (at 9 ^ 90°, i.e. f Ri u —s/2), where the Born asymmetry is maximal, this asymmetry is 
extremely small: 


A° = ^—^ Ri 9.4968-10”®. (8) 

s^ 

It is the main aim of this paper to estimate the contribution of some classes of two-loop contributions, which have 
some logarithmical enhancement. As for the non-enhanced ones - they have an order of {—tjm%){ajTr)'^ rs 10“^^ 
for the CR of MOLLER. Below we consider contribution to the vertex function AV/ for on-mass-shell electrons 
pI = pf = m? and the space-like 4-momentum of the virtual photon = —q^ = — (pi — Pi)'^ ^ rn? in two-loop level 
from the class of Feynman diagrams containing the intermediate states with FF- and Z-bosons (see Fig. [1]). Due to 
vertex renormalization condition AV/|q 2 ^q = 0 the corresponding contribution is proportional to Thus we 

restrict ourselves by the condition 


= i = W,Z. 


(9) 
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III. VERTEX SUBGRAPHS WITH EXTRA W AND Z BOSONS 

The one-loop level expressions for Wm 7 -verte x (see subgraph of process e(pi) —>■ e(pi — k) + j(k) where electron 
with momentum (pi — k) is off-mass-shell in Fig. 1(a)) has a form: 


= -ieu{pi - 

J J \ m^y m^y-k^xxj 


( 10 ) 

( 11 ) 


0 0 


where u}± = 1 ± 75 , g = ejSw and e is the electron charge value (e = |e| > 0). Here and below we use the common 
notation x = l — x, y = l — y, etc. Analogously, one-loop vertex with one additional Z-boson (see corresponding 
subgraph in Fig. |l(b)[ ) looks like: 

^2 


Vu. (pi> k) = -ieu{pi - k)j^{a ± 75 ) u{pi) 


9 


I^{k^) = I dy I dx [ —2 In 


0 0 

Both functions /“’^(fc^) are normalized as J“’^(0) = 0. 

Integrating over 4-momenta of photon k we get the contribution of (flUl) into full two-loop vertex (see Fig. |l(a)[ ) 
as 


— k^xx 


(4cw)^167r2 
k^xx 


i\e) 


m^y rn^zy ~ J 


( 12 ) 

(13) 


AU“ = -ie ^ 

(167r2)2V’ 


ja _ ja _i_ ra _ 


(Tk 


V 


kP-{kP‘ — 2pik){kP‘ — 2p'-^k) 


7T,nk^) 


= u{p[)jx{p'i -k + m) 7 ^(pi -k + m)j^uj_u{pi), 


(14) 

(15) 


where two terms 2 ^ in ([T4| correspond to two terms in (fTTI) . Omitting the terms of order 0{p^ ) we write down 


the contribution as 






f7r2 m'^{k^ — 2 pifc)(fc 2 — 2p'^k) 


(16) 


|fc 2 |<g„ 


Using Feynman parameters trick one can integrate over loop momenta k and obtain the result as a sum of ultra-violet 
finite (UV-finite) and ultra-violet divergent (UVD) parts: 




m 


w 


In^ 

62 


— 1 I dx -f UVD-part, 6 ^ = (pix + p'^x)^ = + Q^xx. 


(17) 


Here and below we use the same notations for unrenormalized and renormalized quantities. Thus after renormalization 
of we obtain expression 


0‘^ 1 777.^ T 

= :;;; 2 -cXrf) 7 ,dW_M(pi), '^i = 3 is ^ 5.0406. 


‘‘W 


(18) 


Here and everywhere below the number value corresponds to CR of MOLLER. 
Second term in da can be written as 


i 2 tp = / yydy / dx 


l- 2 x f d^k 


N° 


1 — X J f7r2 (fc2 — ct2)(/c2 — 2pik){k‘^ — 2p[k) ’ 


(19) 


where cr^ = m^y/ixx). Again using standard manipulations we arrive to 


i y 11 

i 2 fi = u{Pi)itJ^^-u{pi) J yydy J dx^Y^ J J 2yidyi 


0 0 




( 20 ) 
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where Da = Viba + 6^ = {xipi + xip[)'^ = w? + xiXiQ"^ and A is the UV-regularization parameter. Applying 

the renormalization procedure we obtain 


ra — 
^2)1 — 


m 


^C2M(pi)7;,a;_u(pi) 

w 

11 1 


C 2 =- j dxi j 2yidyi j yydy J dx^ 


l-2x 


■ X 


— X 


0 


0 0 0 

, 1 XiXiylxx\ pwxx{l - yiXi){l - yiXi) 

In 1 + 


pwyyi + yfxxxixi 


= -0.0930. 


pw yyi 

The final expression for contribution of bK-vertex subgraph to the vertex function (see Fig. |l(a)[ ) is 

(16^2^2 (^1 + 

Contribution of Z-vertex subgraph (see Fig. |l(b)[ ) has a form 

g’^Sna ^ f d'^k N^, 




k'^{k'^ — 2pik){k'^ — 2p[k) 


AV = -ie- 

(4cu/)2(167r2)2 
= “bi)7A(pi -k + m)7^(pi -k + m)'yx{a ± 75)^w(pi). 

In the similar way we obtain for contribution of Z-vertex subgraph to the vertex function 

AV;? = -ie^ (ci + C2 )u(p'i)7^u(pi), 




m\ (4cu/)^(167r2)2 


where 


2 , 

9“‘'Q2 
1 




4 = -/(Ici J 2yidyi j yydy J dx-^-^x 

0 0 0 0 

/ I xixiylxx\ pzxx{l - yiXi){l - yiXi) 

pz In IH-z- I + 


pz yyi 


pzyyi +yfxxxixi 


= -0.0944. 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


IV. ELECTROWEAK ELECTRON MASS OPERATOR INSERTION TO THE VERTEX FUNCTION 


Now let’s consider the set of Feynman diagrams of vertex type containing electron Mass Operator (MO) with internal 
Z or IF bosons insertions (see Fig. l(c)|l(d)'|. The relevant contribution to the vertex function has a form 


ApMO = -*e|- [F/ + V^] , F; = y 


d'^k 


m 


where i = Z,W and the numerator V* is: 

^ (J, 


with 


A:2(fc2 _ 2p[k) ’ 

Nl = u{p[)jx{p[ - fc + to)7^(pi - fc + TO)7^cW'u(pi)M®(fc,pi), 


(z) ^ 2g^(a±75)^ ^ g'^uj- 

(4cw)^87r2 ’ IOtt^ 

Mass operator M^{k,pi) of electron with both external legs off-mass-shell looks like: 

1 1 


M*(fc,pi) = J xixidxi J 


dz 


ml — Xiz{k'^ — 2pik) 


(26) 


(27) 


(28) 


(29) 
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Standard Feynman procedure of joining the denominators and the loop momentum integrating gives us 

1 1 


y* = / xidxi I ^ I" ^ ^ 


N„, 


z J — 2p[k){k^ — 2pik — af) ’ 


rrii 


= u(p'i)7a(p'i - fc + m)jf^{pi -k + m)'yxu{pi), a‘ = —- 


Xiz 


which can be simplified to the form: 


(30) 


1 111 

= u{p[)jfj,u{pi)V\ ^ / xidxi J ^ J dx J 2ydy (in 


A 2 Q 2 

D, Di 


0 0 0 

After renormalization and expansion on powers of Q'^/ml we obtain 

1 111 

dx / 2ydy 


D, = by + xya^. (31) 


Vi = [ xidxi 

TO- J 

0 0 0 0 


t/' 


- / - 2 - \ 
xizyx + —[y - y xx) 


Final expression for the contribution to the vertex function (see Fig. 1(c) |l(d) ) is 


= ie— 


Q" .z 9^ 


(1 ± a)'^u{p[)jf,u{pi) + Tfz2^W)7i«‘^-'“(Pi) 

w 


““^271 Lto'|"3(4c^) 24^2V—w ^.3 


with 




w _ 11 ///w I 2 _ 2 9925 
3 " 6 g2 + 3 - 


(32) 


(33) 


(34) 


V. CONTRIBUTION OF DIAGRAMS CONTAINING IFIF 7 , IFIF 77 VERTIGES 

Below we consider diagrams containing VFVF 77 , WW^ vertices only because their contributions are associated with 
logarithmic enhancement. Let’s consider the Feynman diagram with virtual photon which is emitted from the initial 
electron and absorbed by the final electron (see Fig. 1(e)). The relevant contribution to the vertex function is 


^ 2(167r2)2>’ 


= 


d'^k 


1 


d'^ki 




iir'^ (/c2 — A2)(fc2 — 2pik){k‘^ — 2p[k) J iir'^ kf{{k + ki — pi)^ — m'^){{k + ki — p()2 — m^) ’ 

where A is the photon mass and 

VcTfiTj — 9(7tJ>i.‘^Pl Pi ^ ^ 1 ) 7 ^ “t“ 9Pi ^ ki'jfj Qrjai^ Pi Pi “t“ 2{k -\- kiYj^j 

Nr,a- = u(p'i)7a(p'i - k)-fnki'Ya{pi - k)j^uj-u{pi). 


(35) 


(36) 


Doing the similar treatment as it was done above one can integrate over loop momentum ki , renormalize the amplitude 
of this subgraph and obtain 


yf = 


3g2 f d*k u{p[)-fx{p'i - fc)7M(Pi “ k)-f^u!-u{pi) 


^ 2m^ J (fc2 — A2 )(/c 2 — 2piA:)(/c2 — 2p(fc) 


(37) 




After integration over k one gets: 


..e 3Q^ 

^ 2 m^ 


dx / 2ydy 


In 


-+ y xx) 


De = yH^ + X^y. 


(38) 


0 0 
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Further we use simple integrals 


= -ln^, 

_|_ X2y Jj2 X'^ ’ 


2ydy 


Q^xx 




dx = 1, 


Q\ , 1 2 Q" 

—TT in —^ ax = in —^ — 


62 


3 ’ 


and obtain 


3Q2 _ 


AK®=*e„ 9 
^ 2m^ 




P = 2 


VK 

2 


47ra^^ ^ 
2(167r2)2 
2 ^2 


In 


”i_ 21 .^,l„ 2 j_ 2 )-ln»«' "'-O'- ” 


<3" 


A 2 


= —2 In ■ 


2 77^2 

In ^ - 40.388. 


A 2 


(39) 


The diagrams in Figs. l(f)|l(g)|l(h) has a general enhancement factor which is associated with the collinear photon 
emission in vertex. Let’s demonstrate this in general. The common structure for all three diagrams contains the 
emission of photon with momentum k from initial electron. This leads to the following structure of the amplitude: 


yf,g,h ^ ^ 

^ 167r2 


(Pk u{p[)0^{pi -k + m)jxu{pi) 
77r2 fc 2 (fc 2 — 2 pifc) 


(40) 


where corresponds to the remaining part of Feynman diagram and different for each diagram. We note that the 
dominant contribution to this integral comes from the integration region of small photon momentum (i.e. |fc^| <C ai^) 
and thus we can omit k in the remaining part of vertex amplitude, containing the momenta of a bF-boson. Joining 
the denominators we have 


yf,g,h _ ^ 

^ 167r2 


dx2xu{pi)pi\ O 




[d^k 
i(pi) / — 


1 


{{k — pix)"^ — m 2 x 2 ) 




and this approximately equals 




i(pi), 


R : 


g\k~xpi 

The diagram containing the JFJF 77 vertex (see Fig. |l(f)[) gives 


L = In 


m 


w 


n2 ■ 


AVf-2ieR^^S 


Pi 


d'^ki 


N- 


Act 


(41) 


(42) 


Z7r2 kj {kf - 2piki - rri^) (fc^ - 2p(fci - m^) ’ 

= u{p'i)lcjhlx^-u{pi), S^^Xa = ‘2gp.„gxa - g^xg^a - g^aPi^X- 
The loop momentum integral do not have ultraviolet as well as infrared divergences. Standard manipulations lead to 


For the Feynman diagram with WW^ vertex shown in Fig. |l(g)| we have 


AV;? = -2iR 


47rag2 f d'^ki u{p[)^nkiJxOJ-u{pi)V-^^’"^V2iJ''^ piv 

327r2 J J7r2 y2 (j^2 _ 2p^ki - {kl - 2p'^ki - m^) 


(43) 


(44) 


where vertices have the form 


= (2pi - k.fg''^^ + {2k, - p,Yg-^ + (-pi - k,Y g^'', 

V 2 Y'" = {-Pi - p'l + 2fci)^p'''^ + (2pi - p'l - kiYgZ + (2p'i - Pi - kiYgY 


(45) 
























Retaining in the nnmerator the terms quadratic over loop momenta one gets 

1 1 


AV^f = j xdx J y^dy^ (~t) ’ 


Finally we have for the contribution of Feynman diagram shown in Fig. |l(g)[ 

<^2 agM3_ , 

L -T:^—u[p^)-/f,u}-u{pi), 


Ays = -ie _ ^ _ 

^ 327r^ 36 


and the diagram shown in Fig. |l(h)| gives the similar result 

ag'^ 67 


Ay = -ie ^ 

^ 647r3 36 


u{p'ihf,uj-u{pi). 


Dg sa ym^. 


(46) 


(47) 


(48) 


VI. NUMERICAL CONTRIBUTION TO THE LEFT-RIGHT ASYMMETRY 


Collecting the result of considered two-loops contributions one can put the total result in the form 

+ AV^'^° -f AV)f+/+9+'‘ = + C^K^, 

where 


= *e-^(l ±a)^ 


ag^ 


(4cvi/)^2 567r^ 


u{jp'i)lg,u(pi), 


K'f = le 


.Q‘^ag'^_ , 


and the coefficients look like 

= -8(c 5 -f 4) + 32cf = 70.5285, 


= -4(c? -f 4) + 8cf -f 3/"’^“ = -312.382. 


(49) 

(50) 

(51) 


(52) 


Let’s note by index C the contributions investigated here, i.e. C = a^b ,..., h. As specific corrections to observable 
parity-violating asymmetry induced by contribution C we choose the contribution to the asymmetry (AA)c and the 
relative corrections D'^-. 


(AA)c = 


\M-c 


-\Ml 


<A|2 


c (AA)c \Mc | 2 -|M ++++|2 
— 


AO 


|A4o-P - |A4(f-+++i2’ 


(53) 


SI-AIqi |yvig I — |yvig 

The physical effect of radiative effects from contribution C to observable asymmetry is determined by the relative 
correction (see [l^ for more details): 


g, _ A<=' - AO _ - 5^ 

- ^5 - i + jc ’ 


(54) 


where the relative correction to unpolarized cross section is 5'^ = o-qq/o-qq. For two-loop effects (where is small) 
the approximate equation for relative correction to asymmetry takes place: 

Contributions to asymmetry of Z and W types are 


(AA)z = -16aC 


zQ^ 


ag^TT 


m\ (4cw)^(167r2)2 

which give the relevant numerical values: 

{AA)z = -2.5410-10■^^ 


(AA)w = 4(7 


^ ag^TT 

(167r2)2 


iAA)w = -3.1983-10 


-10 


(55) 


(56) 


Taking into account that in CR of MOLLER the Born asymmetry A° = 94.97 ppb the numbers for relative corrections 


are 


= -0.0000267, dJ = -0.0033677. (57) 

We can see that effects have the same negative sign, first is rather small, but the second one is at the edge of region 
of planned one per cent experimental error for MOLLER and thus will be important for future analysis of MOLLER 
experimental results. 
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